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Abstract. Radio polarimetric observations of Active Galactic Nuclei can reveal the magnetic 
field structure in the parsec-scale jets of these sources. We have observed the 7-ray blazar 
3C 273 as part of our multi-frequency survey with the Very Long Baseline Array to study 
Faraday rotation in a large sample of jets. Our observations re-confirm the transverse rotation 
measure gradient in 3C 273. For the first time the gradient is seen to cross zero which is 
further indication for a helical magnetic field and spine-sheath structure in the jet. We believe 
the difference to previous epochs is due to a different part of the jet being illuminated in our 
observations. 



1. Introduction 

According to current knowledge, jets of AGN are magnetically launched from a rotating black 
hole or accretion disk [HE1E]- It is then natural to expect the magnetic field structure in the jet 
to be helical [3]. We do not know, however, if the helical magnetic field structure dominates down 
to parsec-scale jets or if the field becomes predominately tangled due to re-collimation shocks 
or interaction with the external medium |5j. Magnetic field structure of parsec-scale AGN jets 
can be studied by using polarimetric multifrequency observations. Evidence for helical magnetic 
field structure could be a detection of a Faraday rotation measure (RM) gradient transverse to 
the jet [6]. 

Faraday rotation occurs when polarized waves propagate through non-relativistic plasma 
within or external to the source [7J. It can be described with a linear dependence between 
the observed electric vector position angle (EVPA, Xobs) an d wavelength squared (A 2 ) by the 



2 Present address: Cahill Laboratory of Astronomy and Astrophysics, California Institute of Technology, 1200 E 
California Blvd, Pasadena CA 91125 



following formula 

e 3 A 2 f 

Xobs = Xo + 87r 2 6Qm 2 c 3 J UeB ' dl = Xo + RMA 2 , (1) 

where xo is the intrinsic EVPA and RM is the rotation measure, related to the electron density 
n e and the magnetic field component B parallel to the line of sight. The sign of the RM then 
depends on the direction of the line of sight component of the magnetic field with RM being 
positive if the magnetic field points towards the observer and negative if it points away from the 
observer. Thus a helical magnetic field viewed directly from the side would have a RM gradient 
going from positive to negative and crossing zero in the middle of the jet. Due to relativistic 
aberration effects blazars that are observed at small viewing angles are viewed from the side in 
the rest frame of the jet. Recent simulations [8] show that launching the jet from a rotating 
black hole can indeed produce such structure. 

The first observational evidence for a RM gradient across a jet was reported in quasar 3C 273 
by Asada et al. [9J, who detected a RM change from 450 to 250 rad/m 2 across the jet in 
observations between 4.6 and 8.6 GHz in 1995. The gradient was confirmed by Zavala & Taylor 
[TO] who in their two observing epochs in 2000 detected a change from 2000 to rad/m 2 across 
the jet of 3C 273 about 5 mas from the core. They attribute the large difference in the RM values 
in comparison to the previous observations [9] to the better resolution of their high frequency 
observations between 12.1 and 22.2 GHz. Follow-up observations on the gradient at frequencies 
between 4.6 and 8.6 GHz in 2002 reveal time variations on time scales of several years fTT] . 
Since these observations were done at the same angular resolution as the first detection [9j, 
the difference between the epochs cannot be due to observational effects. The total rotation 
they detect is larger than can be caused by internal Faraday rotation, as it would cause severe 
depolarization. They also rule out external clouds as the cause of the rotation as the variability 
time scales are too short to cause variations in foreground clouds. They conclude that the 
rotation is most likely caused by a sheath around the jet. RM gradients have been detected in 
other sources as well but the issue remains controversial due to the difficulty in assessing the 
reliability of a RM gradient [12J. Recently, RM gradients have been detected in the kpc-scale 
jets as well, see e.g. Gabuzda et al. (these proceedings). 

In 2006 we conducted a large survey of Faraday rotation in 191 sources in the MOJAVE 
(Monitoring of Jets in Active galactic nuclei with VLBA Experiments) sample (Hovatta et al. 
2011 in preparation). In this paper we report our observations of the RM gradient in 3C 273. 

2. Observations and data reduction 

We observed 3C 273 on March 9, 2006 and June 15, 2006 with the Very Long Baseline Array 
(VLBA) at 8.1, 8.4, 12.1 and 15.3 GHz. The initial data reduction and calibration were 
done in the standard manner as described in the MPS cookboolE The imaging and self 
calibration were done using the Difmap package [13J. For more details see Ref. [Hj. The 
absolute EVPA calibration of the 15 GHz observations were done previously as part of the 
MOJAVE project using the stability of the D-term phases of some antennas and IFs [15l [T6j . 
The 8 and 12 GHz bands were calibrated using nearby observations from the VLA/VLBA 
polarization calibration databas^], and single-dish observations from University of Michigan 
Radio Astronomy Observatory. We estimate the uncertainty in the EVPA calibration to be 3°, 
2° and 4° in the 15, 12 and 8 GHz bands, respectively. 

The (u,v)-range of all the bands was clipped to correspond to each other, and the images at all 
the frequencies were restored to the 8.1 GHz beam. Because the absolute position of the source 
is lost during our calibration process, we used 2D cross-correlation to properly align the images 
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at different frequencies. The images of the different bands were then combined to create the RM 
maps by performing a linear least-squares-fit of the EVPA vs. A 2 in each pixel independently. 
Goodness of the fit was estimated by a x 2 criterion. Details of the procedure will be reported 
in Hovatta et al. 2011 (in preparation). 

3. RM maps of 3C 273 

The RM maps for 3C 273 are shown in Fig. [TJ The top panels show the RM map with color 
indicating the RM in rad/m 2 overlaid on top of total intensity contours at 15 GHz. The second 
panel from the top shows the error in RM, calculated from the variance-covariance matrix of the 
linear fit. The third panel from the top shows the EVPAs at 15 GHz which have been corrected 
for the Faraday rotation. In the bottom panel the sign of the RM is shown. 

The maps are blanked for pixels with absolute RMs of over 20 000 rad/m 2 . These are typically 
concentrated on the core or edge regions of the jets and they do not coincide between the two 
epochs. Therefore we do not think these extreme RM values are real but are artifacts that occur 
by chance due to random noise in the data and utt ambiguities present in the analysis. We 
are currently investigating this in detail with simulations (Hovatta et al. 2011, in preparation), 
and our preliminary results show that it is possible to get these large RM values purely due to 
random noise. The errors in the RM are typically around 70 rad/m 2 rising to as high as 200 
rad/m 2 in the jet edges where the signal-to-noise ratio in total intensity is typically smaller. The 
Faraday-corrected EVPAs are mostly perpendicular to the optically thin jet indicating that the 
magnetic field is predominantly parallel to the jet. The jet region about 7 mas from the core is an 
exception with clearly differing EVPA direction. This could be caused by a shock compressing 
the magnetic field which would also explain the enhanced total intensity and polarized emission 
from this region. The sign map of the RM shows clearly bilateral structure indicative of a helical 
magnetic field as shown by simulations [8j. 

4. Transverse RM gradient 

Our observations in Figs. [2] and [3] show a clear gradient transverse to the jet and along the 
whole jet length. The maximum gradient is seen about 3-7 mas from the core where the RM 
changes from 500 to -600 rad/m 2 . Moving down the jet the gradient becomes less pronounced 
as was also noted by Attridge et al.[17j. Our observations were done at resolution which is in 
between those by Refs. and Ref. [10] and that can explain some of the discrepancy in our 

values. The most notable difference is that our RM gradient goes from positive to negative RM 
values, as expected for a helical field viewed from the side, in contrast to Refs. P HH [TO] where 
the values are always positive. Our positive values also do not extend to as high values as in 
the maps of Ref. [TO] . The most likely explanation is that we are seeing a different part of the 
jet being illuminated by components moving down the jet as our observations are several years 
apart from earlier observations [18] . The effect of seeing a different part of the jet at different 
epochs is clearly demonstrated in the stacked images of the MOJAVE program [14] which show 
how the jet at any given epoch may fill only a portion of the region visible in the stacked image. 
In fact, if we calculate the inner jet direction in 2000 and compare it to the jet direction in 2006, 
the difference is 10°. We are thus seeing a portion of the jet which is 10° to the South of the 
observations in Ref. [10J. This together with the slightly different resolution of our observations 
could explain the difference in the RM values. A similar effect is seen in the RM structure of 
the radio galaxy 3C 120 over several years time scales [T9l 120] . 

Even though our observations are only three months apart we can still see differences in the 
RM between the two epochs. In particular, the core regions strongly differ. This is probably due 
to the complex sub-structure of the core region in 3C 273 [13 [21] . If we ignore the unreliable 
jet edges, the jet RM values are very similar in the two epochs and differences are likely due to 
observational effects such as signal to noise in the differing jet regions resulting in bad A 2 -fits in 
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Figure 1. Left: 
Epoch 2006-03-09. 
Right: Epoch 2006- 
06-15. From top 
down: RM in color 
overlaid on 15 GHz 
total intensity con- 
tours, error of the 
RM, RM-corrected 
EVPA at 15 GHz 
overlaid on the 15 
GHz polarized flux 
contours, sign of 
the RM. In the top 
panel the map has 
been blanked for 
very high absolute 
RM values. 
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one epoch or the other to be in different parts of the jet. The gradient looks very similar over 
the two epochs. The consistent RM gradient over the length of the whole jet suggests a scenario 
involving a sheath around the jet. Recent simulations show that observing a bilateral structure 
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Figure 2. Left: RM map of 3C 273 in 2006-03-09 with additional high-RM pixels of over 4000 
rad/m 2 blanked to show the gradient. Right: RM fit of the core and a jet component 7 mas 
from the core as marked by the dots in the RM map. Bottom: Transverse slices of the jet as 
marked by the lines in the RM map. The absolute EVPA calibration error has been subtracted 
from the error bars because it is the same in every pixel. The beam size along the slice is shown 
as a line in the bottom left corner. 

in the sign map of the RM is a strong indication of helical magnetic field surrounding the jet 
which cannot be easily explained by a random foreground screen [8]. 

Detailed multi-epoch observations of the RM and polarization structure are needed to 
fully understand the Faraday-screen properties, as was recently done for 3C 120 [20]. Their 
observations showed that even though they detect a gradient and variations in the jet RM, these 
are not connected with the polarization changes. This suggests that the Faraday medium may 
be a distant external screen not connected with the jet. 

5. Conclusions 

We have observed 3C 273 with the VLBA as part of our large survey to study the Faraday 
rotation in AGN jets. Our observations confirm the rotation measure gradient transverse to the 
jet. For the first time the gradient is seen to cross zero which is further evidence for a helical 
magnetic field and a sheath around the jet. 
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Figure 3. Same as Fig. [2] but for epoch 2006-06-15. 
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